Criteria for fitting MEG dipoles with fMRI position constraints 
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1 Introduction 

Since MEG and EEG inverse problems are "ill- 
posed” (i.e., solutions cannot be uniquely 
determined), additional information such as fMRI 
and PET activation data can be used to impose 
additional constraints [l]-[6]. For example, fMRI 
data were used to fix the locations of equivalent 
current dipoles at activation areas [1][2][5] or to 
give higher weights for finding dipoles at activation 
areas [3][4][6]. The latter method is capable of 
obtaining dipoles not only at fMRI activation areas 
but also at areas where fMRI fails to show 
activation. It means, however, that the non¬ 
uniqueness nature of the inverse solutions remains in 
this method, and that flexibility in finding dipoles 
and uniqueness of solutions are a tradeoff to be 
optimized. 

In comparison, the former method of fixing dipole 
locations shows active dipoles only at fMRI 
activation areas, but has the advantage that the 
solutions are uniquely soluble by a least-squares 
method with no additional assumptions. This 
method has been applied to the analysis of actual 
data [5]. However, there are two issues that need to 
be clarified further: (1) how to place discrete dipoles 
for approximating spatially extended neural sources, 
and (2) how to account for influences on 
neighboring dipoles, or "crosstalk" when plural 
dipoles are fixed in a small neighborhood. Although 
average behaviors were previously analyzed using 
Monte Carlo simulations [3], their spatially 
anisotropic characteristics need to be investigated 
for effective application of the method. 

In the present work, we carried out simulations of 
solving for dipoles with fixed positions (constrained 
dipoles) to clarify the above two issues, and we 
applied the obtained criteria to actual data. 

2 Methods 

2.1 Simulation 

We used ASA Software Version 2.2 (ANT Software 
BV) to simulate the spatial dependence of fitting 
constrained dipoles (Simulation A), and of crosstalk 
between neighboring constrained dipoles 


(Simulation B). For both simulations, we assumed a 
spherical headshape with a radius of 75 mm, and the 
sensor layout of a 148-channel whole-head MEG 
system (Magnes 2500 WH, 4-D Neuroimaging). We 
defined a coordinate system such that the y-axis 
pointed from the right preauricular point to the left 
preauricular point, the x-axis passed through the 
nasion and was perpendicular to the y-axis, and the 
z-axis pointed upward and was perpendicular to the 
x-y plane. Since radial current produces no magnetic 
fields, we modeled neural sources as dipoles placed 
in a tangential direction (y-direction). We simulated 
source depths of 10, 30, and 50 mm, and source 
extensions including point, linear or planar 
distributions in 20-mm steps inside the headshape. 
We calculated the magnetic fields produced by the 
source dipoles (forward calculation) and added 
random noise to make a signal-to-noise ratio of 3. 

In Simulation A, we fitted magnitudes and 
orientations of a single constrained dipole to the 
source magnetic fields using a least-squares method 
(inverse calculation). We repeated this by varying 
the fitted dipole’s location in x-, y-, and z-directions 
around the source dipole, and evaluated the position 
tolerance which is allowable to obtain a “good fit,” 
defined as a correlation between source and fitted 
dipole magnetic fields exceeding 85%. 

In Simulation B, we fitted two constrained dipoles, 
dipole 1 and dipole 2, to the source magnetic fields. 
Dipole 1 was placed in 4 different positions 10 mm 
from the source in the x, y, +z and -z directions (x 
and y being symmetric about 0). The location of 
dipole 2 was varied within three planes (x-y, y-z, 
and z-x) intersecting each dipole 1 position. To limit 
computation time, dipole 2 positions were simulated 
only along the x-, y-, and z-axes and along 45- 
degree diagonals within the three planes. We defined 
crosstalk to be the ratio of magnitudes of dipole 2 to 
dipole 1, and determined the separation tolerance for 
obtaining a crosstalk of less than 50%. Simulation B 
was carried out for a limited set of source 
parameters: a point source at depths of 10, 30, and 
50 mm, and a 40 x 40-mm extended planar source at 
a depth of 30 mm. 



2.2 Application of the method to actual data 

Based on the above simulations, we obtained two 
criteria (see “Results”) which we then applied to a 
set of actual fMRI and MEG data taken from one 
subject. The experiments were previously carried 
out using a Sternberg paradigm. Subjects were 
prompted by a visual cue to do a phonological task, 
a subvocal rehearsal of 6-character sequences that 
were encoded in short-term memory [7]. 

In the fMRI experiments, neural activity in the 
rehearsal condition was contrasted with that in the 
control condition in which the subjects encoded 
nothing and did no rehearsal. Echo-planar images 
(EPIs) were taken using a 1.5-T fMRI system 
(Magnetom Vision, Siemens AG) with the 
parameters TR = 12.65 s, TE = 66 ms, flip angle = 
90 degrees, voxel size = 2.2 x 2.2 x ( 7 + gap 2.8 ) 
mm. Using MEDx version 3.2 (Sensor Systems, 
Inc.), we motion corrected and detrended the EPIs, 
coregistered them to T1-structural images with a 
voxel size of 4 x 4 x 4 mm, applied a spatial filter 
(FWHM = 8x8x15 mm), did statistical tests using 
cross-correlation methods, and applied a 
significance threshold to individual voxels 
(probability corrected for multiple comparisons of 
10' 2 and uncorrected probability of 6 x 10" 6 ). 

MEG data were taken with the 148-channel whole- 
head system described previously using the same 
stimuli as those used in fMRI experiments. Unlike 
the fMRI experiments, we measured MEG in the 
rehearsal condition but did not subtract neural 
activity in the control condition from that in the 
rehearsal condition. The MEG data were sampled at 
a frequency of 678 Hz, averaged over 150 epochs, 
and filtered with a passband of 1-40 Hz. 

3 Results 

Simulation A showed that position tolerance 
strongly depended on direction. It was narrowest, 20 
mm, for the x-direction which was perpendicular to 
the orientation of source dipoles on the tangential 
(x-y) plane (Fig. 1). The tolerance for the point 
source was similar to that for the extended sources, 
but the latter was a few mm narrower than the 
former. 

Simulation B showed that separation tolerance 
expressed as the coordinate difference between 
dipole 2 and dipole 1 was largest, 37 mm, along an 
axis tilted by -45 degrees in the x-y tangential plane 
((x2-xl) = -(y2-yl), Fig. 2), except that there was no 
tolerance in a radial direction. For simplicity, we can 
treat the maximum separation tolerance to be 
37 mm, because two independent sources separated 


more than 37 mm apart are rarely placed along the 
radial direction. The difference in tolerance between 
the point and the extended sources was less than 1 
mm. Among the three source depths, the largest 
tolerance was obtained for the source depth of 30 
mm. 
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Figure 1: Position tolerance for fixing the location 
of a single constrained dipole to obtain a 
correlation > 85% between the source and 

constrained-dipole magnetic fields (Simulation A). 
The constrained dipole's field was determined after 
fitting the dipole magnitude and orientation to 
simulated source magnetic fields. The constrained- 
dipole location was varied along x-, y-, and z- 
directions around the simulated sources. The graph 
shows the minimum tolerance values for the various 
source distributions. 


„ 40 
E 

£, 20 


0 


0 
o 
c 
CO 

® - 20 | 
o 

H -40 




ill 


x2 - xl 
-H my2-yl 
nz2-zl 
] 0 (x2 - xi)= 
-(Y2-Y1) 


xl yl zl zl 
10 10 10-10 [nnnn] 
Shifts of dipole 1 location 


Figure 2: Separation tolerance for two neighboring 
constrained dipoles which maintains crosstalk at 
less than 50% (Simulation B). The tolerances are 
expressed in terms of the difference in coordinates 
between constrained dipole 2 (x2, y2, z2) and 
constrained dipole 1 (xl, yl, zl) when dipole 1 was 
shifted from the source by 10 or -10 mm in the x-, y- 
and z-directions. The graphs show the results 
obtained for the source depth of 30 mm. 










































From both simulation results, we obtained the 
following two criteria: 

(1) Divide fMRI activation areas if they extended 
more than 20 mm. 

We placed dipoles at each of the fMRI activation 
areas, divided according to (1), and solved for the 
constrained dipoles. Then we applied the second 
criterion: 



Latency [ms] 


(2) Sum up moments of neighboring dipoles if they 
are closer than 37 mm and show simultaneous 
activity. 

We applied these criteria to actual data. By dividing 
the extended fMRI activation areas (Fig. 3), 46 
dipole locations were determined. Since fMRI 
results were obtained by subtraction of activation 
from the control condition (i.e., neural activity of 
doing nothing when prompted by the visual cue), we 
added 6 more activation areas which were assumed 
from other experiments at the calcarine fissures, 
extrastriate visual cortices, and posterior inferior 
temporal areas in both hemispheres. Using ASA 
Software, the magnitudes and orientations of the 52 
constrained dipoles were fitted to the MEG data 
(Fig. 4). We found that the data during a pre-trigger 
period conformed closely to the Rayleigh 
distribution, and that 24 dipoles among the 52 
showed significantly large magnitudes (significance 
level of 1%) after correction for multiple 
comparisons. After summing up the neighboring 
dipoles according to criterion 2, 11 dipole groups 
survived (Fig. 5). They included activity at the left 
posterior superior temporal area (203-210 ms) and 
the left inferior frontal area near Broca's area (247- 
272 and 416-427 ms). These areas have been 
considered to be related to the relevant phonological 
processes [7]. 



Figure 3: One slice of an fMRI activation map which 
contrasted neural activity during subvocal rehearsal 
of 6-character strings and during the control 
condition. 


Figure 4: MEG waveforms averaged over 150 
epochs during subvocal rehearsal of 6-character 
strings. 
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Figure 5: Locations and latencies of 11 dipole 
groups that were obtained using our two criteria. 
Abbreviations: posterior inferior temporal area 
(post.inf.T), intra-parietal sulcus (IPS), posterior 
superior temporal area (post.sup.T), prefrontal area 
(preF, BA 10/11/46/47/45), middle temporal area 
(midT), inferior frontal area (infF), medial inferior 
occipital area (med.inf.Occ). 


4 Discussion 

From the present simulations, two criteria were 
obtained: (1) a minimum position tolerance of 20 
mm for obtaining good dipole fits (correlation > 
85%) and (2) a maximum separation tolerance of 37 
mm for obtaining low crosstalk (< 50%). Their 
strong dependence on direction reflects the fact that 
the magnetic field changes more when dipole 
locations are shifted in a tangential direction 
(especially along the x-direction, i.e., perpendicular 
to the source dipole) than in a radial direction. 

In our approach to fMRI-constrained dipoles, we 
tried to reduce the number of dipoles by dividing the 
extended activation areas into fractions that are as 
large as the position tolerance. This was intended to 
avoid unnecessarily increasing the degrees of 
freedom, which tends to make inverse problems 
underdetermined. Nevertheless, we still needed to 
account for the phenomenon that neighboring 





dipoles placed within the separation tolerance could 
suffer from excessive crosstalk. Since this tolerance 
depended on directions whose exact values were not 
known, we needed to check coincidence of 
activation in latency, in addition to the separation 
between the dipoles, to judge whether we should 
sum neighboring moments or not. 

Applying the two criteria to the set of fMRI and 
MEG data measured using a phonological task, we 
obtained 11 dipole groups with significantly large 
summed magnitudes. Two of the significant dipole 
groups appeared in areas known to be involved in 
phonological processes (L. post.sup.T and L. infF, 
including Broca's area) [7]. Although we applied the 
present methods to only one set of data taken from 
one subject, the results demonstrate that our criteria 
are promising in the analysis of data which includes 
many neural sources. 
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